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Abstract 

Identical neutral kaon pair correlations are measured in y/s = 7 TeV pp collisions in the ALICE 
experiment. One-dimensional K^K^ correlation functions in terms of the invariant momentum dif- 
ference of kaon pairs are formed in two multiplicity and two transverse momentum ranges. The 
femtoscopic parameters for the radius and correlation strength of the kaon source are extracted. The 
fit includes quantum statistics and final-state interactions of the ao/fo resonance. K®K® correlations 
show an increase in radius for increasing multiplicity and a slight decrease in radius for increasing 
transverse mass, tot, as seen in %n correlations in the pp system and in heavy-ion collisions. Trans- 
verse mass scaling is observed between the K^K^ and n% radii. Also, the first observation is made 
of the decay of the f 9 (1525) meson into the K^K^ 1 channel in pp collisions. 



*See Appendix fAlfor the list of collaboration members 
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1 Introduction 

In this paper we present results from a K°K° femtoscopy study by the ALICE experiment [fl] |2l in 
pp collisions at y/s = 7 TeV from the CERN LHC. Identical boson femtoscopy, especially identical 
charged %% femtoscopy, has been used extensively over the years to study experimentally the space- 
time geometry of the collision region in high-energy particle and heavy-ion collisions [3]. Recently, the 
ALICE collaboration has carried out a charged %% femtoscopic study for pp collisions at yfs = 7 TeV. 
This study shows a transverse momentum dependence of the source radius developing with increasing 
particle multiplicity similar to the one observed in heavy-ion collisions H. The main motivations to carry 
out the present K''K° femtoscopic study to complement this %% study are 1) to extend the transverse pair 
momentum range of the charged %% studies which typically cuts off at about 0.8 GeV/c due to reaching 
the limit of particle identification, whereas K°'s can easily be identified to 2 GeV/c and beyond, 2) since 
K° is uncharged, K°K° pairs close in phase space are not suppressed by a final-state Coulomb repulsion 
as is the case of charged %% pairs, 3) K^K^ pairs close in phase space are additionally enhanced by 
the strong final-state interaction due to the ao/fo resonance giving a more pronounced signal, and 4) 
one can, in principle, obtain complementary information about the collision interaction region by using 
different types of mesons. The physics advantage of items 1) and 4) is to study transverse mass scaling, 
of the source size which is considered a signature of collective behaviour in heavy-ion collisions J3J. 
If bit scaling is present, the source sizes from different species should fall on the same curve vs. mj. 
Thus, comparing results from %% and K°K*', particularly at very different mj, would be a good test of 
this scaling. Item 3) can be used as an advantage since the final-state interaction of K°K° via the ao/fo 
resonance can be calculated with a reasonable degree of precision. Previous K°K*' studies have been 
carried out in LEP e + e~ collisions @[6j|7], HERA ep collisions [8], and RHIC Au-Au collisions |9). 
Due to statistics limitations, a single set of femtoscopic source parameters, i.e. radius, R, and correlation 
strength, A, was extracted in each of these studies. The present study is the first femtoscopic K°K^ study 
to be carried out a) in pp collisions and b) in more than one multiplicity and transverse pair momentum, 
ki, range, where kj = \pri + Pti\ /2 and pn and pji are the two transverse momenta of the K^K*' pair. 

2 Description of Experiment and Data Selection 

The data analyzed for this work were taken by the ALICE experiment during the 2010 yfs = 7 TeV pp 
run at the CERN LHC. About 3 x 10 8 minimum bias events were analyzed. 

K° identification and momentum determination were performed with particle tracking in the ALICE 
Time Projection Chamber (TPC) and ALICE Inner Tracking System (ITS) fl]|2]]. The TPC was used to 
record charged-particle tracks as they left ionization trails in the Ne-C02 gas. The ionization drifts up 
to 2.5 m from the central electrode to the end caps to be measured on 159 padrows, which are grouped 
into 18 sectors; the position at which the track crossed the padrow was determined with resolutions of 2 
mm and 3 mm in the drift and transverse directions, respectively. The ITS was used also for tracking. 
It consists of six silicon layers, two innermost Silicon Pixel Detector (SPD) layers, two Silicon Drift 
Detector (SDD) layers, and two outer Silicon Strip Detector (SSD) layers, which provide up to six space 
points for each track. The tracks used in this analysis were reconstructed using the information from 
both the TPC and the ITS ; such tracks were also used to reconstruct the primary vertex of the collision. 
For details of this procedure and its efficiency see Ref . |f2] . 

The forward scintillator detectors, VZERO, are placed along the beam line at +3 m and —0.9 m from 
the nominal interaction point. They cover a region 2.8 < T] < 5.1 and —3.7 < tj < —1.7, respectively. 
They were used in the minimum-bias trigger and their timing signal was used to reject the beam-gas 
and beam-halo collisions. The minimum-bias trigger required a signal in either of the two VZERO 
counters or one of the two inner layers of the SPD. Within this sample, events were selected based on the 
measured charged-particle multiplicity within the pseudorapidity range \rj \ < 1.2. Events were required 
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Fig. 1: Invariant mass distributions of n + n~ pairs in the four multiplicity-^ ranges used in the study. Kj used in 
this analysis were identified by the cut 0.490 GeV/c 2 < m < 0.504 GeV/c 2 . Also shown is a Gaussian+linear fit to 
the data points. 

to have a primary vertex within 1 mm of the beam line and 10 cm of the centre of the 5 m long TPC. This 
provides almost uniform acceptance for particles with | T] | < 1 for all events in the sample. It decreases 
for 1.0 < |t]| < 1.2. In addition, we require events to have at least one charged particle reconstructed 
within \ ri \ < 1.2. 

The decay channel K° — > n + 7i~ was used for particle identification, with a typical momentum resolu- 
tion of ~ 1% ifTOl . The distance of closest approach (DC A) of the candidate K° decay daughters was 
required to be < 0.1 cm. Figure Q] shows invariant mass distributions of candidate K° vertices for the 
four multiplicity-^ ranges used in this study (see below) along with a Gaussian + linear fit to the data. 
The invariant mass at the peaks was found to be 0.497 GeV/c 2 , which is within 1 MeV/c 2 of the accepted 
mass of the K° ifTTl . The average peak width was a = 3.72 MeV/c 2 demonstrating the good K° mo- 
mentum resolution obtained in the ALICE tracking detectors. A vertex was identified with a K*' if the 
invariant mass of the candidate tz + tz" pair associated with it fell in the range 0.490-0.504 GeV/c 2 . As 
seen in Figure [T] the ratio of the K° signal to signal+background, S/(S + B), in each of the four ranges 
is determined to be 0.90 or greater. The minimum K*' flight distance from the primary vertex was 0.5 
cm. Minimum bias events with two of more K°'s were selected for use in the analysis, with 19% having 
greater than two K°'s. A cut was imposed to prevent K°K° pairs from sharing the same decay daughter. 

3 Results 

Figure |2] shows a K°K*' correlation function in the invariant momentum difference variable Q mv = 
yj 1 Q 2 — Qq, where Q and Qo are the 3-momentum and energy differences between the two particles 
respectively, for all event multiplicities and kj. The correlation function was formed from the ratio of 
"real" K ( $ 'K° pairs from the same event to "background" K°K° pairs constructed by event mixing of ten 
adjacent events. Bins in Q- mv were taken to be 20 MeV/c which is greater than the average resolution of 
<2inv resulting from the experimental momentum resolution. Also, the enhancement region in <2mv of the 
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Fig. 2: Inclusive (all event multiplicities and k^) K^K^ Qinv correlation function. Plotted in the insert to the figure 



is the invariant K^K^ mass distribution, dAVdm(KjKj), in the vicinity of the small peak at Qi m ~ 1.15 GeV/c. 



correlation functions for source sizes of ~ 1 fm is ~ 200 MeV/c. Thus the smearing of the correlation 
function by the experimental momentum resolution has a negligible effect on the present measurements. 
The three main features seen in this correlation function are 1) a well-defined enhancement region for 
Ginv < 0.3 GeV/c, 2) a non-flat baseline for Q mw > 0.3 GeV/c and 3) a small peak at Q- mv 1.15 GeV/c. 

Considering feature 3) first, fitting a quadratic + Breit-Wigner function to the invariant K°K° mass dis- 
tribution, d./V/dm(K^K°), around this peak, where m(K''K°)= 2yJ (<2inv/2) 2 +nt^, we obtain a mass of 

1518 ± 1 ± 20 MeV/c 2 and full width (T) of 67 ± 9 ± 10 MeV/c 2 (giving the statistical and systematic 
errors, respectively). This is plotted in the insert to Figure [2] Comparing with the Particle Data Group 
meson table IfTTl . this peak is a good candidate for the f^(1525) meson (m = 1525 ± 5 MeV/c 2 , T = 73+!j 
MeV/c 2 ). This is the first observation of the decay of this meson into the K°K° channel in pp collisions. 

In order to disentangle the non-flat baseline from the low-<2i nv femtoscopic enhancement, the Monte 
Carlo event generator PYTHIA fT2l[T3l was used to model the baseline. PYTHIA contains neither quan- 
tum statistics nor the K ( S 'K° — > ao/fo channel, but does contain other kinematic effects which could lead 
to baseline correlations such as mini -jets and momentum and energy conservation effects [4|. PYTHIA 
events were reconstructed and run through the same analysis method as used for the corresponding ex- 
perimental data runs to simulate the same conditions as the experimental data analysis. The PYTHIA 
version of the invariant mass distributions shown for experiment in Figure Q] yielded similar S/(S + B) 
values. As a test, the K°K° background obtained from event mixing using PYTHIA events was compared 
with that from experiment. Since the background pairs do not have femtoscopic effects, these should ide- 
ally be in close agreement. A sample plot of the experimental to PYTHIA ratio of the background vs. 
<2inv is shown in Figure [3]for the range A^ c h 1-11, &t < 0.85 GeV/c. The average of the ratio is normalized 
to unity. It is found that PYTHIA agrees with the experimental backgrounds within 10%. The Monte 
Carlo event generator PHOJET lfT4l [TBI was also studied for use in modeling the baseline, but it was 
found to not agree as well with experiment as PYTHIA, and thus was not used. 
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Fig. 3: Ratio of K^K|? experimental background to PYTHIA background vs. Q{, 
GeV/c. The average of the ratio is normalized to unity. 



for the range N ch 1 - 1 1 , k T < 0.85 



K°K^ correlation functions in Q mw were formed from the data in four ranges: two event multiplicity 
(1-11, > 11) ranges times two kj (< 0.85, > 0.85 GeV/c) ranges. Event multiplicity was defined as the 
number of charged particles falling into the pseudorapidty range |t]| < 0.8 and transverse momentum 
range 0.12 < pj < 10 GeV/c. The two event multiplicity ranges used, 1-11 and > 11, correspond to 
mean charged particle densities, (dVch/drj), of 2.8 and 11.1, respectively, with uncertainties of ~ 10%. 
PYTHIA events were used to estimate (dN c h/dri ) from the mean charged-particle multiplicity in each 
range as was done for Table I of Reference [4], which presents ALICE Tin results for pp collisions at 
y/s = 7 TeV, event multiplicity having been determined in the same way there as in the present work. This 
is convenient since the K°K° source parameters from the present work are compared with those from the 
%% measurement below. About 3 x 10 8 minimum bias events were analyzed yielding about 6 x 10 6 K°K° 
pairs. A similar number of PYTHIA minimum bias events used for the baseline determination were 
also analyzed. This was found to give sufficient statistics for the PYTHIA correlation functions such 
that the impact of these statistical uncertainties on the measurement of the source parameters was small 
compared with the systematic uncertainties present in the measurement. The femtoscopic variables 7? and 
X were extracted in each range by fitting the experimental correlation function divided by the PYTHIA 
correlation function with the Lednicky parametrization [9 ] based on the model by R. Lednicky and V.L. 
Lyuboshitz lfl6l . This model takes into account both quantum statistics and strong final-state interactions 
from the ao/fo resonance which occur between the K°K° pair. The K° spacial distribution is assumed to 
be Gaussian with a width R in the parametrization and so its influence on the correlation function is from 
both the quantum statistics and the strong final-state interaction. This is the same parametrization as was 
used by the RHIC STAR collaboration to extract R and A from their K''K° study of Au-Au collisions 
0. The correlation function is 



C(anv)=AC , (2inv) + (l-A) 
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Fig. 4: Experimental and PYTHIA K®K® correlation functions for the four multiplicity-^x ranges. 



f(k*) is the s-wave K°K° scattering amplitude whose main contributions are the s-wave isoscalar and 
isovector fo and ao resonances (9], R is the radius parameter and A is the correlation strength parameter 
(in the ideal case of pure quantum statistics A = 1). a is the fraction of K°K° pairs that come from the 
K°K° system which is set to 0.5 assuming symmetry in K° and K° production [9|. As seen in Eq. (2), 
the first term is the usual Gaussian from quantum statistics and the second term describes the final-state 
resonance scattering and both are sensitive to the radius parameter, R, giving enhanced sensitivity to this 
parameter. The scattering amplitude, f(k*), depends on the resonance masses and decay couplings which 
have been extracted in various experiments The uncertainties in these are found to have only a small 
effect on the extraction of R and A in the present study. An overall normalization parameter multiplying 
Eq. (1) is also fit to the experimental correlation function. 

Figure|4]shows the experimental and PYTHIA K ( $ 'K° correlation functions for each of the four multiplicity- 
kj ranges used. Whereas the experimental correlation functions show an enhancement for Q mv < 0.3 
GeV/c 2 , the PYTHIA correlation functions do not show a similar enhancement. This is what would 
be expected if the experimental correlation functions contain femtoscopic correlations since PYTHIA 
does not contain these. PYTHIA is seen to describe the experimental baseline rather well in the region 
Ginv > 0.4 GeV/c 2 where it is expected that effects of femtoscopic correlations are insignificant. Fig- 
ure [5] shows the experimental correlation functions divided by the PYTHIA correlation functions from 
Figure |4] along with the fits with the Lednicky parametrization from Eqs. (l)-(3). The fits are seen to 
qualitatively describe the correlation functions within the error bars, which are statistical. 

Figures [6] and [7] and Table [TJ present the results of this study for A and 7? parameters extracted by fitting 
the Lednicky parametrization to K°K° correlation functions as shown in Figure [5] The source parameters 

are plotted versus mj = J (&t) 2 + m K to observe whether mj scaling is present (as discussed earlier) and 
statistical + systematic error bars are shown. The statistical uncertainties include both the experimental 
and the PYTHIA statistical uncertainties used to form the correlation functions, as shown in Figure |4] 
The largest contributions to the systematic uncertainties are 1) the non-statistical uncertainty in using 
PYTHIA to determine the baseline and 2) the effect of varying the <2inv fit range by ±10%. These 
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Fig. 5: Experimental K^K, correlation functions divided by PYTHIA correlation functions for the four 
multiplicity-^T ranges with femtoscopic fits using the Lednicky parametrization. 



were found to be on the order of or greater than the size of the statistical uncertainties, as can be seen 
in Table [T] The method used to estimate the systematic uncertainty of using PYTHIA was to set the 
PYTHIA KvK° background distribution equal to the experimental background distribution in the ratio of 
correlation functions, e.g. forcing the ratio plotted in Figure[3]to be exactly unity for all Q my . The ratio of 
correlation functions then becomes the ratio of the experimental to PYTHIA real pair distributions, which 
is then fit with the Lednicky parametrization to extract the source parameters. Parameters extracted from 
these correlation functions were then averaged with those from Figure [5] and are given in Figures |6]and 
[7] and Table Q] This method is similar to that used in estimating systematic uncertainties in other K°K° 
measurements JUIH. 

To see the effect of the tfo/Zo final-state interaction (FSI) term in the Lednicky parametrization, the 
correlation functions in Figure [5] were fit with Eqs. (l)-(3) for two cases: 1) quantum statistics + FSI 
terms, i.e. a = 0.5 in Eq. (2), and 2) quantum statistics term only, i.e. a = in Eq. (2). Case 2) 
corresponds to the usual Gaussian parametrization for R and A. The results of these fits are shown in 
Table (2) Including the FSI term in the fit is seen to significantly reduce both R and A, i.e. R by ~ 30% 
and A by ~ 50%. The FSI is thus seen to enhance the correlation function for Q mv — > making A appear 
larger and making the enhancement region narrower resulting in an apparent larger 7?. A reduction in R 
and A when including the FSI term was also observed, but to a lesser extent, in the STAR Au-Au K°K° 
study [9]. A larger effect of the ao/fo resonance on the correlation function in pp collisions compared 
with Au-Au collisions is expected since the two kaons are produced in closer proximity to each other in 
pp collisions, enhancing the probability for final-state interactions. 

Within the uncertainties, the mj dependence of A is seen in Figure [6] to be mostly flat with A lying at an 
average level of ~ 0.5 — 0.6, similar to that found in the ALICE Tin results for pp collisions at y/s = 7 
TeV [4]. In Tin studies the A smaller than 1 has been shown at least in part to be due to the presence 
of long-lived meson resonances which distort the shape of the source so that the Gaussian assumption, 
which the fitting functions are based on, is less valid ifTTl . This same explanation is possible for the 
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Fig. 6: A parameters extracted by fitting the Lednicky parametrization to K^K^ correlation functions as shown in 
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Table 1: K^K^ source parameters from Lednicky fits for y/s = 7 TeV pp collisions. Statistical and systematic 
errors are shown. 



kj range 


N ch range 


<*r> 


(cW ch /dT]) 


A 


R 


(GeV/c) 




(GeV/c) 






(fm) 


<0.85 


1-11 


0.52 


2.8 


0.66 ±0.07 ±0.04 


0.99 ±0.04 ±0.04 


>0.85 


1-11 


1.32 


2.8 


0.48 ±0.04 ±0.06 


0.75 ±0.02 ±0.07 


<0.85 


> 11 


0.52 


11.1 


0.53 ±0.07 ±0.05 


1.15 ±0.05 ±0.05 


>0.85 


> 11 


1.32 


11.1 


0.54 ±0.04 ±0.09 


1.00 ±0.02 ±0.07 



Table 2: K^K^ source parameters comparing a — 0.5 (quantum statistics+FSI) and a = (quantum statistics only) 
fits to Figure fusing Eqs. (l)-(3). Statistical errors are shown. 

kj range N c h range A R A R 

(GeV/c) (fm) (fm) 

a = 0.5 a = 0.5 a = a = 

<0.85 1-11 0.64 ±0.07 0.96 ±0.04 1.36±0.15 1.35 ±0.07 

>0.85 1-11 0.50± 0.04 0.81 ±0.02 1.07 ±0.09 1.05 ±0.04 

<0.85 >11 0.51±0.07 1.12 ±0.05 0.97±0.15 1.64±0.11 

>0.85 >11 0.56±0.05 1.03 ±0.02 0.89±0.10 1.37 ±0.07 



present A parameters extracted from the K''k[' correlation functions. For example, the and K* mesons 
with full widths of T ~ 4 and T ~ 50 MeV/c 2 , respectively, could act as long-lived resonances compared 
with the extracted source scale of R ~ 1 fm, the larger scales being unresolved in the first few <2inv bins 
but still depressing the overall correlation function. 

In Figure |7] the dependence of the extracted radius parameters on the transverse mass and event multi- 
plicity are shown. Also shown for comparison are R parameters extracted in the same event multiplicity 
ranges from a tctc femtoscopic study by ALICE [4] in 7 TeV pp collisions. Looking at the mj dependence 
first, the K°K° results alone suggest a tendency for R to decrease with increasing mj for both multiplicity 
ranges. The %% measurements also show this decreasing trend for the high multiplicity range, but show 
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Fig. 7: R parameters extracted by fitting the Lednicky paramettization to K®K® correlation functions as shown 
in Figure [5] Also shown for comparison are R parameters extracted in the same event multiplicity ranges from a 
%% femtoscopic study by ALICE |4| in pp collisions at ^fs = 7 TeV. Statistical (darker lines) and total errors are 
shown. The highest mj pion N c h > 1 1 point and lower mj kaon N c h > 1 1 point have been shifted by 0.03 GeV/c 2 
for clarity. 



the opposite trend for the low multiplicity range, R increasing slightly for increasing mj. Taken with 
the %% results the K°K*' results for R extend the covered range of mj to ~ 1.3 GeV/c, which is more 
than twice the range as for 7171. The lower mj points for K°K° which are in close proximity in mj to the 
highest mj points for %% are seen to overlap within errors, showing mj scaling. The mj dependence of 
R combining both particle species is seen to be weak or nonexistent within the error bars. Looking at the 
multiplicity dependence, a tendency for R to increase overall for increasing event multiplicity is seen for 
both %% and K*'K'' measurements as is observed in %% heavy-ion collision studies ifTBl . 

The multiplicity-/Mx dependence of the pion femtoscopic radii in heavy-ion collisions is interpreted 
as a signature for collective hydrodynamic matter behaviour Q. The corresponding measurements in 
pp collisions at ^Js = 7 TeV show similar behaviour [4J. However, important differences with heavy- 
ion collisions remain, for example the low-multiplicity dependence of R in pp for pions seems to 
increase with increasing rather than decreasing as with heavy-ion collisions, as already mentioned 
earlier. The interpretation of these pp results for pions is still not clear, although model calculations exist 
that attempt to explain them via a collective phase created in high-multiplicity pp collisions |[T9ir20ll2lTl . 
If such a collective phase is hydrodynamic-like, the mj dependence of the radii should extend to heavier 
particles such as the K*' as well, as shown in Ref. ETTl . The measurements presented in this paper provide 
a crucial cross-check of the collectivity hypothesis. The interpretation is, however, complicated by the 
fact that in such small systems particles coming from the decay of strong resonances play a significant 
role l22l ; simple chemical model calculations show that this influence should be relatively smaller for 
kaons than for pions. So far, no model calculations are known in the literature for any KK correlations 
in pp collisions for mj > 0.7 GeV/c 2 , but the results measured in the present study should act as a 
motivation for such calculations. 



4 Summary 

In summary, identical neutral kaon pair correlations have been measured in y/s = 7 TeV pp collisions in 
the ALICE experiment. One-dimensional K°K° correlation functions in terms of the invariant momen- 
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turn difference of kaon pairs were formed in two multiplicity and two transverse momentum ranges. The 
femtoscopic kaon source parameters R and X have been extracted. The fit includes quantum statistics and 
final-state interactions of the ao/fo resonance. K°K^ correlations show an increase in R for increasing 
multiplicity and a slight decrease in R for increasing mi as seen in %% correlations in the pp system and 
in heavy-ion collisions. The universality of the mj dependence of the extracted radii, i.e. mj scaling, is 
also observed within uncertainties for the K°K"' and nn radii. 
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